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Abstract 

We study the impact of the QCD DGLAP evolution on the geometric scahng 
of the gluon distributions which is expected to hold at small x within the satu- 
ration models. To this aim we solve the DGLAP evolution equations with the 
initial conditions provided along the critical line = Ql{x) with Ql{x) ~ 
and satisfying geometric scaling. Both fixed and running coupling cases are stud- 
ied. We show that in the fixed coupling case the geometric scaling at low x is 
stable against the DGLAP evolution for sufficiently large values of the parameter 
A and in the double logarithmic approximation of the DGLAP evolution this hap- 
pens for A > ANcas/i^- In the running coupling case geometric scaling is found 
to be approximately preserved at very small x. The residual geometric scaling 
violation in this case can be approximately factored out and the corresponding 
form-factor controlling this violation is found. 



1 Introduction 



Perturbative QCD predicts very strong power-law rise of the gluon density xg{x, Q'^) 
in the hmit a; — > 0, where, as usual, x denotes the momentum fraction carried by the 
gluon and is the scale at which the distribution is probed. This strong rise can even- 
tually violate unitarity and so it has to be tamed by screening effects. Those screening 
effects are provided by multiple parton interactions which lead to the non-linear terms 
in the (BFKL and/or DGLAP) equations pi] - |]13|. These non-linear terms reduce the 



growth of gluon distributions and generate instead the parton saturation at sufficiently 
small values of x and/or ^ 



[|0 



Increase of the gluon distribution and emergence of the saturation effects imply 
similar properties of the measurable quantities which are driven by the gluon, like 
the deep inelastic structure function F2{x,Q'^). This can be most clearly seen in the 
dipole picture of deep inelastic scattering in which the virtual photon - proton total 
cross section a^*p{x, Q^) ( a^*p{x, Q"^) ~ -^2(2;, Q'^)/Q^)) is linked with the cross section 
adp{x,r) describing the interaction of the qq colour dipole with the proton, where r 
denotes the transverse size of the dipole M, |21|, E2l E5|. The dipole-proton cross section 



is determined by the gluon distribution in the proton and in leading order approxima- 
tion we just have adp{x,r) ~ as{l/r'^)r'^xg{x, 1/r^). Increase and/or saturation of the 
gluon distribution in the small x limit implies similar increase and/or saturation of the 
dipole-proton cross section and of the cross-section cT^.p(x, Q^). 



The successful description of all inclusive and diffractive deep inelastic data at 
HERA by the saturation model [| 



( see also and [^) suggests that the screening 
effects might become important in the energy regime probed by the present collid- 
ers. Important property of the dipole cross section which holds in this model is its 



a. 



geometric scaling, i.e. dependence upon single variable r = r'^Ql{x) where Qs{x) 
is the saturation scale. This leads to the geometric scaling of a-y*p{x,Q'^) itself, i.e. 
x.Q"^) = f {Q"^ / Ql{.x)) which is well supported by the experimental data from 
Geometric scaling of the dipole cross-section should imply similar scaling 
of the quantity as{Q'^)xg{x, This type of scaling is also found to be an intrinsic 

[E0||. It turns out that for 



HERA 26 



property of the non-linear evolution equations 
the equations of type 



k) 
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UgK - as4>^{x, k) , ( 



Nra, 



a. = 



IT 



(1) 



where is a linear evolution kernel (for example of BFKL type) there exist a region 
in X and k space such that 

(f){x,k) = (j){Ql{x)/k^) for k^ < Ql{x) . (2) 

For example in the case of the Balitsky-Kovchegov |12| equation, where K is the 
BFKL kernel, the saturation scale Ql{x) has been found to have a general power like 



2 



dependence on x , Ql{x) = Q^x'^. The coefficient A which is approximately equal 
to in this case, is then a universal quantity and does not depend on the initial 
conditions for the evolution ||T6| - [pO| . 

The main purpose of this paper is to analyse possible compatibility of this scaling 
with the DGLAP evolution equations. It is expected that the non-linear shadowing 
effects should be weak in the region 'to the right' of the critical line defined by the satu- 
ration scale Q^{x) , i.e. for Q"^ > Q^{x), see Fig. |l|. In order to study possible impact of 
the DGLAP evolution we shall therefore assume the geometric scaling parametrisation 
along the critical line and inspect the structure of the solution of the DGLAP equation 
with those initial conditions. This way of providing the initial conditions along the 
critical line = Q^{x) rather than at = Qq with x independent reference scale Qq 
is the characteristic feature of the saturation effects 0]. 

The content of our paper is as follows. In the next section we give semianalytical 
insight into the solution of the DGLAP equation with the starting distributions pro- 
vided along the critical line. We study separately the fixed and running coupling cases. 
In Section 3 we present numerical analysis of our solutions and finally in Section 4 we 
give our conclusions. 

2 Solution of the DGLAP equations from the start- 
ing distributions provided along the critical line 

We wish to understand possible effects of the DGLAP evolution on the geometric scal- 
ing at low X. This scaling means that certain quantities controlling deep inelastic 
scattering at low x, like the dipole-proton cross section adp{x,r = 1/Q) or the virtual 
photon-proton cross section cr^*p, which are in principle functions of two variables, de- 
pend upon the single variable Q/Qs{x). The saturation scale Qs{x), which also specifies 
the critical line increases with decreasing x 



Let us assume that: 

1. For < Q'^{x) the linear evolution is strongly perturbed by the nonlinear effects 
which generate geometric scaling for the dipole cross section adp{x,r = 1/Q) and 
for the related quantities. 

2. Geometric scaling for the dipole cross-section implies geometric scaling for 
C(s{Q'^)xg{x, where g{x, Q"^) denotes the gluon distribution. This follows 
from the LO relation between the dipole cross section and the gluon distribution, 
i.e. (t(x, r^) ~ r'^as{)-/r'^)xg{x, 1/r^). 

3. Geometric scaling for as{Q^)xg{x, Q"^) /Q"^ holds at the boundary = Q^{x). 




(3) 
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4. For > Q'^{x) the non-linear screening effects can be neglected and evolution 
of parton densities is governed by the DGLAP equations. 

We wish to study possible effects of the DGLAP evolution upon the geometric scaling 
in the region > Q'^{x) after solving the linear DGLAP evolution equations starting 
from the gluon distribution satisfying this scaling and defined along the critical line 
(see point 3 above). We shall discuss the fixed and running couphng cases 
separately. 



2.1 The fixed coupling case 

Let us consider standard leading order evolution of the gluon distribution xg{x, Q"^) 

dxg{x, Q"^) as dz 



^/ -P,,{z)xg{x/z,Q'), (4) 

ZTT Jx Z 



91n(QVA2) 

where, as usual, the Pgg is the gluon-gluon splitting function. For simplicity we have 
neglected possible contribution of the quark distributions. In the moment space this 
equation has the following form 

^ T^lggi^jg^iQ ) , (5) 



ain(QVA2) 2n 
where we have defined the Mellin transform to be 



g^{Q')= ['dxx^g{x,Q^), (6) 

<^ 



and the gluon anomalous dimension is defined as 



7,,(u;) = f'dzz-Pggiz). (7) 



The solution of equation (|) is straightforward and given by 



2\ 1^799 H 



g^(Q^)=g,(uj)(%] . (8) 



We will now seek equation for the moment function go{uj) using the following initial 
condition 

^xgix,Q' = Q'six)) = ^r'x-\ (9) 

where Q^x) is given by equation @. The parameter r° specifies the normalisation 
of gluon distribution along the critical line. This boundary condition follows from the 
geometric scaling condition of the dipole proton cross section a^pir = 1/Q,x) which is 
proportional to asXg{x,Q'^)/Q'^. 
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In order to find solution for go{uj) we use tlie inverse Mellin transform 

xg{x,Q') = 7^ f duox--g^{Q^) , (10) 

where tlie integration contour should be located to the right of the singularities of 
duiQ'^) ill the uj plane. Inserting in equation ( ]T0| ) the DGLAP solution (||) for QuiiQ"^) 
we get 

xg{x,Q^) = —Jdujx'-g,{u:)\^] . (11) 

We now set = Q'i{x) with the saturation scale Ql{x) defined by equation (^, and 
require the geometric scaling initial condition along the critical line = Ql{x) (see 
eq. (|^)). From equations (|^),(0) and (|TT|) we get 

^ J dujgo{uj)x^'^-^^^^^^'^'> = r°x"^ . (12) 

This equation can be regarded as the equation for the function go^uj), i.e. for the 
moment of the gluon distribution at the (x independent) scale Q^. In order to solve 
this equation we take the moment of both sides of equation (|T^), i.e. we integrate both 
sides of this equation over dx for < x < 1 with the weight x'^^~^ and get 

/ duj. — - = . 13 

J [ui-uj- A^7gg(a;)] - A 

We now change the integration variables 

Z=U + \-^lgg{uj) , (14) 
/TT 

which after inversion specifies the function uj{z). Equation ([l3|) in the new variable z 
takes then the following form 

±- f ^^M^qoitMl = ^° . (15) 

27ri J dz {uji — z) lui — X 
We can easily perform the contour integration in Eq.(|T5]) and get 

duj(z) , , , X 

— ^ go{uj{z = UJ,)) = . (16) 

ClZ z=uii LUi — A 

We still need to solve this equation for goioj) and in order to do this we write 

^1 = w + A^733(^) , (17) 
and finally from Eq.(|l^) we obtain 



ga{uj) = [1 + A 



r 







as d-^ggju) 

2vr duj \uj + \^^gg{uj)-\] 
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which defines the solution for go{uj). 



In what follows it is convenient to use directly redefined function go{z) 

Uz) ^ ^goHz)) , (19) 
az 

where from equation (|TB]) we see that 

The solution of the DGLAP equation with the initial condition specified by equation 
(y) then reads 

.3(..0^) = 5-/«w(^] (21) 

where the integration contour is located to the right of the singularities of go{z) and of 
u!{z). If the leading singularity is a pole of go{z) at ^ = A then the leading contribution 
to xg{x, Q"^) at small x is given by 



x^(x,Q^)^r°x-^(^^j , (22) 

where 

uo = cj(A) . (23) 

It should be noted that defines position of the pole of go{uj). In general we have 
ujq < X. From equation ( p2D we get the following leading small x behaviour for the 
gluon density ^xg{x,Q'^)/Q'^ 

a,xg{x,Q') a, ( ^^ fJ^-(-o)-i 

(24) 



2n g2 Ql27i\Ql{x)J 

which respects the geometric scaling i.e. is a function of only one combined variable 
Q'^/Ql{x). Violation of this scaling by the contribution of the (branch point) singularity 
of uj{z) is a non-leading effect at low x. 

The requirement that the pole of g^{z) at 2; = A is the leading singularity imposes 
certain constraints upon A. In general they are difficult to be found exactly since 
the inversion of equation ([I^) cannot be performed analytically when using complete 
form of '-)gg{u)). Analytic solution of equation (|I1|) is however possible in the double 
logarithmic approximation in which 'Jgg{uj) = 'y£,^{uj), where 



2Ar, 



7,T(^) = ' (25) 
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We also have 



is the most singular in ^ part of the gluon anomalous dimension (u;). In this 
approximation we get 

uj{z) = , (26) 

where 

ots = . (27) 

71 

ujo = . (28) 

The condition that the pole of g^{z) at z = A is the leading singularity, i.e. that it 
is located to the right of the branch-point singularity of u^z) at z = 2a/^7A gives the 
following constraint upon the parameter A 

A > 4a, . (29) 

For A < itts the leading singularity is the branch point of uj{z) at z = 2y/dsX and the 
geometric scaling becomes violated. 

It may be interesting to confront our results for the fixed coupling with the proper- 



ties of the exact solution of the non-linear Balitsky - Kovchegov equation |]20|. In this 
case geometric-scaling holds for < Ql{x) and the non-linear effects can be neglected 
for > Q^x). The parameter A specifying the critical line is however not an inde- 
pendent quantity and depends upon the (fixed) coupling a,. In the double logarithmic 



approximation it is given by A = Aag- It follows from equation (p9[) that this is a 
limiting value of the parameter A for the geometric scaling to hold asymptotically in 
the small x limit and so for A = 4a g we expect violation of this scaling for > Ql{x) 
down to the very small values of x [EO]. 



2.2 Running coupling case 

We now pass to the more realistic case with the running coupling. In this case the 
evolution equation for the moment function takes the form: 



2^ 



91n(gVA2) 27r 
where the running coupling in the leading order is given by 



2tt ln(QVA 

with 



2\ ' y^^^) 
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with Nf being number of flavours. In this section we consider only gluonic channel 
therefore we set Nf = 0. The solution of equation (^) reads 



gM) = 9o{uj) { B^iTT^ ) . (33) 

From the above solution we obtain 



2^ 



2. «"« )-^*'-nKiMJ ■ 

and so the result for the gluon distribution xg{x, Q"^) in x space reads in this case 



27r 

where 



(35) 



/o(^) = ^^goiu;) . (36) 
We now impose the geometric scaling condition (|^) onto this solution to get 

which is an equation for foiyj). Solution of this equation is complicated, i.e. exact 
solution generates complicated (branch point) singularity of /o(co') at = A. The only 
observation which we can make is that it should generate behaviour softened by 
inverse powers of the ln(l/x). In order to make some insight into what is going on we 
have to make some approximations. To be precise let us make the approximation by 
setting = A in the argument of iggiyj) that gives 



Aln(l/a;) 



r°x-^ . (3^ 



Making the same approximation in the inverse Mellin transform (RH) we get the 



solution 

27r ^ ^ - ^ " \HQl/K^) ) + ln(gg/A2) ) ' 
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Multiplying and dividing by Q^x' 



we finally obtain 



2ti 



-xg{x, Q^) = X 



HQ^x^/Ql 



^ln(Q2/A2) + Aln(l/x; 



+ 1 



(40) 



The factor proportional to ln(l/x) in the denominator of the expression on the r.h.s. 
of (^) generates violation of the geometric scaling. Thus in the case of running of the 
coupling as{Q^) the scaling behaviour gets violated, it is possible however to factor out 
the effect of this violation. We can also rewrite Eq.(^OD by using the definition of the 
saturation scale and the running coupling to get 



a,iQ')xg{x,Q') Ql 



X 



2ti 



Qi 



1 + 



a,(g2(a;)) 



27r6 



HQ^/Qlix)) 



(41) 



where we see that the violation is proportional to the value of the running coupling 
evaluated at the saturation scale. Consequently when x <^ 1 that is when Qs{x) ^ 1 
the geometric scaling is restored, provided of course that as{Ql{x))\Yi{Q'^ /Ql{x)) ^ 1 
as well. This condition is equivalent to \yi{Q'^ / Ql{x)) <^ \n.{Ql{x) / Is?) . The same 
condition defining the region in which the geometric scaling holds above the saturation 
scale has recently been found in ref. 



27 . 



3 Numerical results 

In this section we present numerical results for the evolution of ordinary DGLAP equa- 
tions for the integrated gluon distribution function with special boundary conditions 
set on the critical line Ql{x) as described in Sec.l. 



3.1 Fixed coupling case 

We start with the simplest case which is the fixed strong coupling. We assume also in 
the first approximation the DLLA limit that is we only keep the most singular part of 
the Pgg splitting function in our simulation i.e. 

27V 

Paai^) = ^^ Nc = 3, (42) 
which results in the following form for the anomalous dimension of Eq. (|^) 

7,,(^) = ^ . (43) 

UJ 

The initial condition for the evolution of the gluon density is assumed to be of the 
form (H). We take A = 0.5 and = 0.1. In Fig.| we show the results of the calcula- 
tion in this case. We illustrate the scaling behaviour of the gluon density by plotting 
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xg{x,Q'^)IQ'^ versus scaling variable r = Q'^/Ql{x) for different values of rapidity 
Y = Inl/x. From Eq.(|2^) we see that this function should scale with r = Q'^/Ql{x). 
The geometric scaling would correspond in this plot (Fig. |^) to the perfect overlap of 
all curves for different values of F , so that they would form one single line. We see that 
up to a good accuracy this function does not depend dramatically on Y and thus on x. 
We do however observe that there is some violation of the scaling at large x. This is 
due to the fact that the geometric scaling expression defined by equation ( P^ is only 
expected to hold asymptotically in the small x limit. At finite x this leading behaviour 
is perturbed by the non-leading contribution given by the branch-point singularity of 
uj{z) at z = 2y/a^ , (see eq. ( p6| ) ). 



To illustrate better the scaling and its violation we have plotted xg{x, Q"^) / versus 
scaling variable r = Q"^ IQl{x) using double-logarithmic scale, see Fig. One clearly 
sees that with increasing rapidity Y the curves do not change and reach asymptotic 
straight line. We have also selected the very low x range of Fig. ^a, which is Fig. 
One can see that in this case the geometric scaling is nearly preserved (we see nearly 
single line for different rapidities). 

The behaviour oixg{x, Q'^) /Q"^ versus r = Q'^ /Ql{x) is clearly governed by a power 
law, with a power which we estimated to be approximately —0.77. From equation (p^ 
and (pSD , and using the values of A and as quoted above we get that the power should 
be ^7gg(i^o) ~ 1 = —0.74 which is in a very good agreement with numerical result. 

Let us note that in the case of DLLA (|^, ujq is a solution of the quadratic equa- 
tion and is given by (0). As previously noticed the real solution exists only for 
A > Xmin = with = asNc/n. We have numerically checked that for A < A^m 
our solution no longer exhibits geometric scaling. It is interesting to note, as we have 
already observed at the end of Sec. 2.1, that exactly the same value of A = A^m for 
a power of saturation scale was obtained from the studies of the nonlinear Balitsky- 
Kovchegov equation [|l^, [1^ performed in [jl6| - pO . 



We next abandon the DLLA approximation and consider more general case with the 
full gluon-gluon splitting function Pgg which gives the following anomalous dimension 



Jggiu) = 2NJ--^— + ^—-^—-^E + ^-H^+2) 

Lc<JCo'+lco' + 2co' + 3 12 



(44) 

where ip is Polygamma function. In this case equation (IT^p with z = X can no longer be 



solved analytically and has to be analysed numerically. However, one can get insight 
into the allowed values of A by making the expansion of the anomalous dimension 
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around u = 0. In this case '^gg{uj) / {2N^ ~ + Ai{Q) + 0{uj) where Ai(0) - 



Using this approximation in ([1^) one finds that now geometric scaling will hold if the 
following condition is satisfied 

4a 

^ - = [l-a-Ji(0)P • ^^^^ 
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We have checked numerically that above approximation works very well and gives very 
close results to the solution of (14) with full uj dependence of anomalous dimension 



In Fig-lja we plot xgi^x.Q"^) /Q"^ as a function of scaling variable Q^/Ql{x) in the 
case of calculation with the full anomalous dimension (@). We have taken A = 0.5, 
and (is = 0.1. We see that the function exhibits geometric scaling (although there is 
some residual violation at larger values of x). The calculated value of exponent from 
numerical calculation is —0.85 which is again in nearly perfect agreement with the 
analytical estimate based on the approximation described above which gives —0.86. 
We also present in Fig. Wo the calculation in the case of A = 0.3 which is below the 



critical value (^5[) equal in this case Xmin = 0.33 for «<, = 0.1. We clearly see that the 
geometric scaling is never present in that case. 

One can study the scaling and its violation in a more quantitative way by examining 
the following expression 

1 dh{Y, r) 



A(r, 



h dY 

where 



fixed 

a. 



(46) 



h{Y,T)^^xg{x,Q') . (47) 

Derivative A{Y, r) should vanish in the region where geometric scaling is satisfied. 
Consequently its deviation from zero will characterise the scaling violation of solution 

We present the quantity A(Y, r) in Fig.§ for the case of calculation with complete 
anomalous dimension and two selected values of A : 0.3 and 0.5. Derivative A{Y,t) in 
Fig.^ therefore illustrates the scaling and its violation for the solution shown in Fig.^. 
From Fig.|^ it is clear that in A = 0.5 the scaling is always reached, even for high values 
of r , that is very far right to the critical line. On the other hand in the case A = 0.3 
the derivative A{Y,t) never vanishes meaning that for A < Xmin the solution of the 
DGLAP equation with the fixed coupling case does not exhibit geometric scaling. 



3.2 Running coupling case 

We consider now the case in which as is running and study the impact of scaling bound- 
ary condition (|^) onto the evolution. We consider full expression for the anomalous di- 
mension in that case 7gg(cu) given by Eq.(P^. The running of the coupling requires that 
the evolution is taken in the region well above the Landau pole. In our case this means 
that one has to evolve with > Ql{x) and we would like to have Ql{x) big enough 
for all values of x. For the purpose of illustration we take that Ql{x) = Ql{x / Xq)~^ 
where = 1-0 GeV^ and Xq = 1.0. This means that at a: = 1 the saturation scale is 
equal = 1.0 GeV^. This assumption might seem artificial considering the present 
phenomenology of lepton-nucleon scattering which suggests that saturation scale could 
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be of the order of 1 GeV^ at a; ~ 10 ^ for the most central colUsions at HERA colUder 



22i , |28[| . However, we use it here for the purpose of illustration of basic effects of the 
evolution with special scaling boundary conditions. We concentrate ourselves here on 
presenting general properties of the solution rather than trying to describe the exper- 
imental data. We also take Nf = 0, that is we are considering pure gluonic channel. 
In Fig.||a we present the results of the calculation by plotting as{Q'^)xg{x, Q'^)/Q'^ ver- 
sus scaling variable r = Q'^/Ql{x) in the case with full gluon anomalous dimension. 
For comparison we also show the calculation performed in the DLLA approximation 
Fig. We see that the geometric scaling is mildly violated in the running coupling 
case, more strongly in the DLLA approximation due to the faster evolution. This fact 



can be understood on the basis of Eq. f^lf ) where the numerical value of the exponent 
of expression on the r.h.s. is much bigger in the DLLA case: bjgg{X = 0.5) = 0.18 in 
the case with full anomalous dimension and h^gg{\ = 0.5) = 13/11 in DLLA case. 

We have tried to estimate whether the violation is consistent with the analytical 
prediction of formula (pOD. In Fig. ^a we present the same quantity as in Fig. ^ but 
multiplied by the scaling variable r = Q'^/Ql{x). The solid black curves in Fig. ^a 
from the upper to lower are for decreasing values of x. One can see that the solution 
exhibits some small violation of the geometric scaling and that the magnitude of this 
violation is smaller for smaller values of x ( the curves are becoming closer and closer 



as X decreases ). This is consistent with general behaviour predicted by equation (pOD 
where the scaling violating factor on the r.h.s. tends to unity when ln(l/x) ^ 1. We 
stress that the observed scaling violation is very small in this kinematical regime. For 
example at very high value of r = 10^ the violation of the scaling is about 5% in a 
huge rapidity range from F = 6 to F = 46. 



It follows from equations ( flOD and (^) that the violation of the geometric scal- 
ing can be approximately factored out. We checked this approximate prediction by 
considering the quantity 



with VF{x 



as{Q^)xg{x,Q^)/Ql{x)VF{x) 
HQ^/Qlix)) 



ln(Q2/A2) + Aln(l/a;; 



+ 1 



(48) 



which according to equation (|^) should be constant with respect to r = Q'^ /Ql{x). 
The results for the above quantity are shown in Fig. |^ (which is Fig. ^a multiplied by 
VF{x) ) where now we see that approximately the geometric scaling is nearly restored 
(curves form a very narrow band) at high values of rapidity . 

Also in the case of running coupling we have studied the features of the geometric 
scaling using the method of the derivative, see Eq. (^B]). The results are shown in 
Fig.^] where it is clear that there is always a region where the geometric scaling is 
(approximately) preserved in the running coupling case, even at very high values of r. 
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This is consistent with formula (|41| ) provided we have asiQl 
with the conclusions of ref 



27, 29 



Inr < 1, 

We have also illustrated in Fig.l 



a; <^ 1 and also 
H the sensitivity 



of the results on the variation of the normalisation for the saturation scale i.e. Qq. 
Changing the parameter Qg from 1 - upper plot in Fig.^ to 0.1 GeV^ - lower plot in 
Fig.^, influences the size of the violation of the scaling. One can see that the geometric 
scaling is postponed to the higher values of rapidity. 



4 Summary and conclusions 



In this paper we studied effects of the DGLAP evolution upon the geometric scal- 
ing. We solved the DGLAP evolution equation for the gluon distribution with the 
initial condition respecting the geometric scaling and provided along the critical line 
= QKx). In the case of the fixed QCD coupling we obtained analytic solution of the 
DGLAP equation with those boundary condition, Eq. (pi]). We also showed that for 
sufficiently large values of the parameter A defining the critical line this solution of the 
DGLAP equation preserves the geometric scaling for the leading term at small x, (see 
Eq. (p2D). In the double logarithmic approximation of the DGLAP equation this hap- 



pens for A > Adis, where as is defined by equation (^). Geometric scaling is however 
violated by effects which are subleading at small values of x. We have also obtained 
approximate solution of the DGLAP equation with the running coupling starting again 
from the boundary conditions respecting geometric scaling along the critical line. In 
the running coupling case geometric scaling is mildly violated for arbitrary values of 
the parameter A yet this violation can be approximately factored out. The size of this 
small violation is controlled by the quantity (^^(Q^) InQ^/Q^. Thus in the region where 
X <C 1 and InQ'^/Ql <C InQ^/A^ the geometric scaling in the running coupling case is 
preserved. Results of the detailed numerical analysis confirmed all those expectations. 
We conclude that the geometric scaling is a very useful regularity following from the 
saturation model. We believe that it might be interesting to incorporate this 'DGLAP 
improved' geometric scaling in the phenomenological analysis of the data. 
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In 1 / X 



Saturation region 
Geometric scaling : 
f(x,Q) = f(Q/q(x) )i 
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critical line Q = Q„ (x) 
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Figure 1: Phase diagram in (In l/x, In Q/A) space. Thick line is a critical line 
= Q^x) which divides the saturation - scaling regime (to the left) and the linear - 
DGLAP regime (to the right). 
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Figure 2: Function xg{x,Q'^)/Q'^ in DLLA fixed coupling case plotted versus scaling 
variable r = Q'^/Ql{x) for different values of rapidities Y — \q.1/x, from Ymin — 6.0 
to Ymax = 46.0 (solid curves from lower to upper) in steps AY — 2 . Dashed curve is 
the input distribution ~ 1/r. 



17 




Figure 3: Function xg{x,Q'^)/Q'^ in DLLA fixed coupling case plotted versus scaling 
variable r = Q^/Ql{x) for different values of rapidities Y — Inl/x. Solid curves - 
solutions, dashed curve - input distribution ~ 1/r. On upper plot (a): solid curves 
from lower to upper are for Y rapidities ranging from Y„iin = 6.0 to Ymax = 46.0 in 
steps AY = 2. Lower plot (b): rapidities ranging from Ymin = 26.0 to Ymax — 46.0 in 
steps AY = 2. 
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Figure 4: Fixed coupling case with complete gluon anomalous dimension •jgg^u). Func- 
tion xg{x,Q^)/Q^ plotted versus scaling variable r = Q^/Ql{x) for different values of 
rapidities Y = \nl/x from Ymin = 6.0 to Ymax = 46.0 in steps AY = 2. Upper plot 
(a): scaling exponent X = 0.5, lower plot (b): scaling exponent X — 0.3 . 
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Y=ln 1/x 



Figure 5: The derivative A(Y, r) from Eg. j^^) as a function of rapidity Y for various 
values of the scaling variable r. Upper plot - scaling exponent A = 0.5, lower plot - 
scaling exponent A = 0.3. 
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Figure 6: The solution as{Q'^)xg{x, Q"^) /Q^{x) in the running coupling case. Rapidity 
range from Ymin = 6.0 to Ymax = 46.0 in steps AY = 2. Upper plot (a): case with full 
anomalous dimension jgg{uj), lower plot: (b): case in DLLA approximation. 
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Figure 7: The solution as{Q'^)xg{x,Q^) /Ql{x) in the running coupling case. We have 
selected high rapidity range from Ymin = 18.0 to Ymax = 46.0 in steps AY = 2. Upper 
plot (a): as{Q'^)xg{x, Q"^) /Ql{x) , lower plot (b): as{Q'^)xg{x,Q'^)/Ql{x)VF{x) where 
the factor VF{x) is a scaling violation factor defined in equation ^4^). 
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Figure 8: The derivative A{Y,t) defined in Eg. as a function of rapidity Y for 
different values of t and different choices of normalisation Q^. The scaling exponent 
A was set to be equal A = 0.5. 
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